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Abstract 


Triangular  factorization,  solution  to  Linear  equations,  inversion,  computation  of  rank  profile 
and  inertia  ( in  the  Hermitian  case)  etc.  of  general  n  x  n  matrices  require  G i  1  operations. 
For  certain  structured  matrices  including  Toeplitz  and  Hankei  matrices  the  computational 
complexity  is  known  to  be  0(n2)  or  better.  These  structured  matrices  often  arise  in  a 
■vide  variety  of  area*  including  Signal  processing,  Systems  theory  and  Communications. 
Fast  li.e.  0(nJ))  algorithms  for  these  structured  matrices  have  been  actively  studied  for 
over  twenty  five  years.  However  almost  all  the  authors  have  assumed  that  the  underlying 
matrices  are  strongly  regular  i.e.  every  principal  submatrix  is  nonsingular.  Althougii  some 
fast  algorithms  have  recently  been  developed  for  certain  problems  involving  some  of  these 
structured  matrices  which  may  have  one  or  more  zero  minors,  several  other  problems  remain 
unresolved;  also  a  unified  approach  towards  these  problems  is  lacking.  In  this  dissertation, 
we  obtain  several  new  results  through  a  unified  approach  to  the  problems  mentioned  earlier. 

First  we  derive  fast  (i.e.  0(n2))  procedures  for  computing  a  ■’modified"  triangular  fac¬ 
torization  which  is  a  LDU  factorization  where  L  is  lower  triangular  (resp.  V  is  upper 
triangular)  with  u.  diagonal  entries  and  D  is  a  block  diagonal  matrix  with  possibly  vary¬ 
ing  block  sizes.  Only  strongly  nonsingular  matrices  will  always  have  a  purely  diagonal 
nonsingular  D  matrix.  For  the  matrices  we  study  the  diagonal  blocks  are  also  structured; 
in  particular  they  are  Hankei  (resp.  Quasi- Toeplitz)  for  Hankei  (resp.  Toeplitz)  and  Quasi- 
Hankei  (resp.  Quasi- Toeplitz)  matrices. 

A  particular  application  of  our  result  is  a  fast  method  of  computing  rank  profile  and 
inertia,  leading  to  alternative  proofs  of  certain  results  due  to  lohvidov  ( 197-1)  on  rank  profile 

and  inertia  of  Hankei  and  Toeplitz  matrices. 

Next  using  the  results  on  modified  triangular  factorization  we  extended  the  Schur  com¬ 
plement  based  approach  of  Chun  [Chu89]  for  inversion  of  strongly  regular  Toeplitz  and 
Hankei  matrices  to  Hermitian  Toeplitz,  Quasi-Toeplitz,  Hankei.  and  Quasi-Hankel  matrices 
with  arbitrary  rank  profile. 

We  apply  our  procedures  for  computing  inertia  to  derive  0(n2)  procedures  for  comput¬ 
ing  root  distribution  of  an  n-th  order  polynomial  with  respect  to  the  unit  circle  and  the 
imaginary  axis.  In  particular  it  has  been  possible  to  derive  the  first  general  recursive  pro¬ 
cedure  for  determining  the  root-distribution  of  a  polynomial  with  respect  to  the  unit  circle: 
the  classical  Schur-Cohn  test  for  this  problem  fails  in  the  presence  of  certain  "singularities." 
which  correspond  to  a  Quasi-Toeplitz  Bezoutian  matrix  being  non-strongly-regular. 


